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An abelian gauge symmetry can be spontaneously broken near a black hole horizon in anti-de 
Sitter space using a condensate of non-abelian gauge fields. There is a second order phase transition 
between Reissner-Nordstrom-anti-de Sitter solutions, which are preferred at high temperatures, and 
symmetry breaking solutions, which are preferred at low temperatures. 



PACS numbers: ll.15.Ex, 04.25.Dm, 11.25.Tq 



I. INTRODUCTION 

I have previously argued that an abelian gauge sym- 
metry can be spontaneously broken by a condensate of 
a charged scalar that forms near the horizon of a non- 
extremal charged black hole, presumably indicating some 
form of superfluidity or superconductivity [l|, Q • Unpub- 
lished numerical calculations based on the model pro- 
posed in indicate that "superconducting" horizons arc 
thermodynamically preferred below some non-zero crit- 
ical temperature. But there are enough parameters in 
the lagrangians considered in P, 0] that it is challenging 
to characterize the degree of universality or robustness 
of such numerical results. I will therefore consider an 
analogous phenomenon in a theory whose lagrangian is 
mostly determined by symmetry principles: 
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2^ 



where 



and 



ge^'^A^Al 



(2) 



(3) 



is the field strength of an SU{2) gauge field. I use mostly 
plus signature, and e'^'"^ is the totally antisymmetric ten- 
sor with e^^^ = 1.1 aim to persuade the Aj^ and gauge 
bosons to leap out of the horizon and condense near it, 
thereby breaking a [/(I) symmetry associated with A'^. 
This should be possible if the gauge coupling g is large 
enough. 

It is convenient to represent the gauge field as a matrix- 
valued one form: 



A ^ Ay^dx" , 



where r° = cr"/2i, so that 



[T ,T \ = e T . 



(4) 
(5) 



(By (T° I mean the usual Pauli matrices.) I will restrict 
attention to the following ansatz: 



'2a 



-hdt^ + dxl + dxl 



dr^ 



A = ^T^dt + w{T^dXl + T^dX2) , 



(6a) 
(6b) 



where a, h, and w arc functions only of r. The gauge 
field (j6bp is a slight simplification of the ansatz considered 
in 0| A substantial literature has grown up around 
similar solutions of Einstein- Yang-Mills theory: see 0, 
for reviews. 



The electrostatic potential $ must vanish at the hori- 
zon for A to be well-defined as a one-form, but I do not 
require it to vanish at infinity: thus the black hole can 
carry charge under the U{1) gauge symmetry generated 



(1) by 



The condensate wij^dxi +T'^dx2) breaks the ^7(1) ro- 
tation symmetry in the X1-X2 plane as well as the C/(l) 
gauge symmetry generated by , but it preserves a com- 
bination of the two. I will require w to be normalizable 
in the sense of making a finite contribution to the norm 



\\A\? 



9 9 ; 



(7) 



where rn is the location of the horizon. Normalizability 
of w is what it will mean for the condensate to form 
"near" the horizon. It is an appropriate requirement in 
the context of studying spontaneous symmetry breaking. 

In summary: the only conserved quantities associated 
with the black hole should be its mass density and its 
electric charge density; and w, if it is non-zero, is a 
condensate whose presence spontaneously breaks U{1) of 
spatial rotations times U{1) of into a diagonal sub- 
group. I want to know when this spontaneous symmetry 
breaking occurs. 
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II. SYMMETRY BREAKING SOLUTIONS 

Plugging ((6a|) and (|6bp into the equations of motion following from ([2]) results in four second order equations of 
motion: 

1 2 <j,2 2 

a" + a'2 + VV + ^e-6-^==0 (8a) 
$" + 2a'$'-2.g2e-4aZ!!_ ^0 (8c) 

together with a zero-energy constraint, 

12e2-a'2 + Ae^-^a' -^ + ^- 2uP - 2g\-^-^-^ + g^e-^'^^ = . (9) 
h L-'h h h'^ h 

I 



TABLE I: Scaling symmetries of the equations of motion (|8]), the constraint ([9]), of series coefficients appearing in (|10p and 
and of the thermodynamic variables appearing in (|13|) . 







dx 


dr 


e'^ 




h 


hi 


Ho 


Hi 




$1 


Po 


Pi 


w 


Wq 


Wi 


L 


9 




/ 


e 


s 


T 


P 




J 


A 





1 





-1 


-1 


2 


2 


2 


-1 











-1 


-1 


-1 


-2 














-3 


-2 


-1 


-2 


-1 


-2 


B 


1 


1 


-1 


-1 


-1 





1 





-3 


-1 





-1 


-2 


-1 


-1 


-2 














-3 


-2 


-1 


-2 


-1 


-2 


C 


1 


1 


1 











-1 











-1 

















1 


-1 


1 





-3 


-2 


-1 


-2 


-1 


-2 


D 
























































1 


-2 


-2 


-2 





-2 





-2 



The metric and the gauge field must be regular at the 
horizon, which 1 assume to be at r = 0. This implies 
series expansions around r = of the form 

a = ao + air + a2r^ + ■ ■ ■ 
h = hir + hiT^ + . . . 

(10) 

$ = <^^r + <i>2r + . . . 

w = Wo -\- wir + W2r'^ + . . . . 

ft is also straightforward to obtain the asymptotic be- 
havior at large r: 

a = log - + ao + • • • 

h = Ho + H^e-^^ + H^er^'^ + . . . ^^^^ 
^ = Po + Pie-°- + . . . 
w = Wie-" + ... . 



tions ([8]) and ([9]) have three scaling symmetries that act 
non-trivially on the parameters. These symmetries are 
summarized in the first three rows of table [H where as- 
signing a charge a to a quantity X means 

A:->A"X. (12) 

Recall that w is required to be normalizablc at infinity. 
This requirement is imposed in the last equation of (jlip . 
but it is not the typical behavior of w: usually there is 
a constant term at large r, or else there is a singularity 
at finite r. The solutions with normalizablc w form a co- 
dimension one locus in the three-dimensional space of pa- 
rameters modulo scaling symmetries. This locus includes 
as one branch the well-known Reissner-Nordstrom-anti- 
de Sitter solutions (hereafter RNAdSV, which have the 
form ^ with w = 0: see for example Q. 



Once the six quantities oq, hi, (f>i, wo, g, and L are 
specified, all the other coefficients in pUj) may be com- 
puted, and one can then use the series expansions to 
specify Cauchy data for the equations of motion ([5]) at 
some radius r slightly outside the horizon. Of these six 
parameters, only three are meaningful, because the equa- 



Thc energy density, entropy density, temperature, 
chemical potential, and charge density, as well as an or- 
der parameter, J, to be discussed further below, can be 
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read off as foflows: 
H3 



2tt 



H = 



P 



K^LHo 
Po 

Pi 



,2ao 



47r 



2ao 



hi 



(13) 



J = 



Energy density, temperature, and chemical potential are 
measured in reference to the Killing time \/Hot rather 
than t itself. The quantities in can be regarded as 
characterizing a thermal state of a dual conformal field 
theory, along the lines of the gauge-string duality d, d, 
1 10( 1 . Their normalizations accord with the conventions of 
[7|. All dependence on k can be removed by defining 
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(27r)3L2 
2 

up 
(27r)3L2 



(27r)3i2 

At2j 

(27r)3L2 



(14) 



where the factors of 27r are included for later convenience. 
Any relation among the quantities in (|13p and (jl4p must 
be expressible in terms of ratios which are invariant un- 
der the scaling symmetries summarized in tabic HI For 
example, the RNAdS solutions have 



e 



1 



7r2p2 



(15) 



When w 7^ 0, one may ask what fraction q of the elec- 
tric charge density is carried by the non-abclian gauge 
bosons outside the horizon. The ratio of the fiux of the 
r3 electric field through the horizon to the flux at infinity 
is 1 — (7. Therefore, 



q = l + Le2""V^— ■ (16) 



q by itself is invariant under the scaling symmetries. 



The conformal field theory dual to the AdS^ back- 
ground under consideration has currents satisfying 
an SU(2) current algebra, where m runs over the t and 
Xi directions. The symmetry breaking that arises from 
non-zero w corresponds to expectation values 



cx J5f , 



(17) 



III. SUMMARY OF RESULTS 

As explained following ((T^), there is a two-parameter 
family of black hole solutions with a non- vanishing, nor- 
malizable condensate, once all possible scaling symme- 
tries are used. A convenient choice of parameters is gL 
and Tj^fp. Other scaling invariants include J//5, q, and 
A//p3/2^ where 



A/ 



(27r)3L2 



(18) 



Here A/ is the difference in the free energy density, e—Ts, 
between a symmetry breaking solution and the RNAdS 
solution with the same T and p. A/ < means that the 
symmetry breaking solution is preferred. 

I studied solutions to ([5]) numerically for about 1600 
different choices of parameters: see figure [TJ I restricted 
attention to solutions with w everywhere positive. Solu- 
tions exist in which w has nodes, but they are probably 
always thermodynamically disfavored because spatial os- 
cillations in w increase energy density. A shooting strat- 
egy was employed to implement the constraint of nor- 
malizability on w. Solutions with non-zero, normalizable 
w were found only for T less than a critical temperature 
Tc- This Tc is, within numerical error, the temperature 
at which the RNAdS solution admits a static linearized 
perturbation, with w non-zero but infinitcsimally small. 
The shooting algorithm was designed to work well near 
Tc- In practice, it worked well for T > Tc/2. 

Tc goes to zero at g = gc ~ 0.8/L and does not ap- 
pear to have singular behavior there; however, the nu- 
merical problem becomes more difficult at small T^, so 
the value of g^. should be regarded as approximate. For 
g ^ 4(7c, scaled thermodynamic quantities such as J / p, q, 
and A///53/2 exhibit nearly universal behavior as func- 
tions of T/Tc from T = down at least to ^T^. This 
universality is related to a large g limit where the back- 
reaction of the gauge field on the metric can be neglected. 

For Q.75Tc <T <Tc,l found good fits to the following 
scaling forms: 



J r 

p 



(19) 



where t = 1 — T/Tc, and qi and ji/2 are positive and 
depend on gL. More approximate fits can be found over 
the same temperature range to the scaling form 



p3/2 



(20) 



where i runs over the values 1,2. The tensor 6f exhibits 
the locking of a spatial [/(I) and a gauge U{1). pT|) 
describes a form of long-range order which infrared fluc- 
tuations probably destroy; however, fluctuations arc sup- 
pressed in the supergravity approximation, where k <C L 
|17j . I will therefore ignore them. 



where /2 is positive and depends on gL. A few solutions 
close to Tc appeared to have A/ very slightly positive. 
But closer examination of some of these solutions using 
highly accurate numerics showed that in fact they have 
A/<0. 

In conclusion: there is a second order phase transition, 
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with simple critical exponents, from Reissner-Nordstrom- 
anti-de Sitter solutions to solutions with normalizable 
non-abelian condensates. The critical temperature scales 
as the square root of charge density and also depends 
on the gauge coupling. The solutions with condensates 
break an abelian gauge symmetry, suggesting some link 
or analogy to superconductivity. An analogy has been 
already suggested between charged anti-de Sitter black 
holes without symmetry breaking condensates and the 
pseudogap state of high Tc materials [ll|, [H, ■ It 
would be desirable to test the validity of this analogy 
more closely, and to extend it if possible. Given the way 
spatial rotations enter into the description of the con- 
densate ([T7)) . it might be possible in some variant of the 
construction I have exhibited to find an analog of p-wave 
or c?-wave superconductivity. Any such analog with high 
Tc superconductivity should be viewed with caution if 



not skepticism, because the underlying degrees of free- 
dom of duals to AdSi vacua are typically large TV gauge 
theories, seemingly distant from semi-realistic construc- 
tions like the Hubbard model. And yet, one might hope 
that dynamics related to the global symmetries (usually 
U (2) or S'0(4), and sometimes larger groups) of the Hub- 
bard model and its relatives can be mimicked, to some 
extent, by non-abelian gauge fields on the gravity side of 
the gauge-string duality. 
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APPENDIX A: STATIC PERTURBATIONS 

The purpose of the appendices is to consider two Umit- 
ing cases of the equations ([8]): the first, in this appendix, 
is the case where w is treated as a small perturbation to 
an RNAdS solution; and the second, in appendixiBl is the 
case where the entire non-abelian gauge field is treated 
as a small perturbation to the Ad5'4-Schwarzschild geom- 
etry. 

The RNAdS background is 



ds^ = -fdt 
f 

$ = pLK 



.^2 dr^ 



1 



4,.2 
1 



r 
L2 



(Al) 



with w = 0. By definition, is the most positive zero 
of /. If If is treated as a static linear perturbation, then 
the only equation of motion we must solve is ((8d|) with 
the w'^ term neglected: 



w H — -w 



f 



p 



0. 



(A2) 



Although (|A2p does not appear to be solvable in terms 
of known functions, it is easy to make a scries expansion 
near the horizon: 



w 



(A3) 



Once wq is specified, all other coefficients in (|A3p can be 
computed. Because (|A3P is linear, it suffices to consider 
the case = 1, even though our real interest is in so- 
lutions with very small w. In addition to (jA3p . there is 
also a solution that has a logarithmic singularity at the 
horizon, but I discard it based on the usual argument 
that the horizon is a non-singular location. Far from the 
horizon, the typical behavior of w is for it to be asymptot- 
ically constant, but by adjusting rjj/L one can arrange 
for w to have 1/r falloff, which is normalizable in the 
sense described around ([7]). The discretely many allowed 
values of rn/L depend on e, p, g, and L, but because of 
scaling symmetries (essentially, the ones in table |l|, di- 
mensionless quantities such as T/ ^/p depend only on gL. 
Thus one may convert the allowed values of th / L into 
allowed values of at some specified gL. The largest 

allowed value of Tj^fp corresponds to a solution where 
w has no zeroes, and in this case the temperature is pre- 
cisely the Tc discussed in the paragraph following ([18]). 
The next-to-largest allowed value of T/ ^fp corresponds 
to a solution with one node; the next-to-next-to-largest 
to two nodes; and so forth. See figure [H 

The last term in (|A3[) shows that the condensate devel- 
ops for essentially the same reason as in the scalar case 
[2|: w acquires a negative effective mass squared near 
the horizon. When the vT" term in (jSdp becomes signifi- 



cant, it decreases or limits the tendency of w to condense. 
An analogous situation in the lagrangian studied in [2l 
would be to have a positive ji/'l^ term in the scalar poten- 
tial, similarly limiting the tendency of the charged scalar 
■0 to condense. 



APPENDIX B: STRONG COUPLING LIMIT 

The large gL limit is simple because the gauge field 
doesn't back-react significantly on the metric [l8|. To 
see this, note that upon scaling A Aj g and F F/g, 
the action for the gauge fields is just -^{F^^Y. The 

stress tensor acquires a similar factor of l/g^, and it is 
this factor that suppresses the back-reaction. Approxi- 
mate solutions to the full equations of motion ([8|) and 
the constraint ([9]) can therefore be constructed by start- 
ing with the j4dS'4-Schwarzschild line element. 



— [-{l-r]j/r^)de^dx^]+-^- 



L2 dP 



r^/r^ ' 
(Bl) 

and then solving the Yang-Mills equations in that back- 
ground. Let us set r/f = 1. Then (jBip has the form ([5a|) 
with 



a = log ■ 



h = l- 



1 



Let us also define 



$ gL^^ 



gL w. 



(B2) 



(B3) 



Then the Yang-Mills equations (|5c|) and (|8dp become 



r{r^ 



1 + 2r'^ 



l>2w - 



r{r^ — 1) 



= 0. 

(B4) 

A charming feature of (jB4p is that there are no free pa- 
rameters in the equations. Factors of g are absent be- 
cause the definitions (jB3|l include the aforementioned 
scaling A A/g. It is perhaps surprising that one also 
needs explicit factors of L in (|B3|) to avoid them in (|B4[) . 
This is because of the choice of radial variable in (|B1[) : if 
I had used z = L'^/r, such factors would not be needed, 
because then L enters into the line element only as an 
overall factor, and overall factors don't effect conformally 
invariant equations like the classical Yang-Mills equations 
of motion. 

The equations (|B4p cannot be solved analytically, but 
they are straightforward to solve numerically, starting 
with expansions 



l>:==$i(r-l) + l'2(r-l) 



w 



Wq + Wi{r — I) + 



(B5) 
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FIG. 2: Left: Temperatures where static linear perturbations of tlie RNAdS solution arise. The curve labeled "no back-reaction" 
is given analytically by (jBSj. It comes from neglecting the back-reaction of "1? as well as of w on the metric. As discussed in 
section |B] neglecting back-reaction is a good approximation when gL is large. Right; Static normalizable perturbations, scaled 
so that m = 1 at the horizon. The solution that is everywhere positive corresponds to the dot at gL = 6 on the T = Tc curve 
in the left-hand plot, and the other solutions correspond to the dots on the other curves. 




near the horizon and 



Starting from (|B3p and (fT4| . it is easy to show that 



Wi 
r 



El 
r 



(B6) 



far from it. The only free parameters are $i and wq: all 
higher coefficients of (|B5p can be determined in terms of 
them. The normahzable form of shown in (jB6[) is not 
the typical behavior: instead, id usually asymptotes to 
a non-zero constant at large r. Requiring normalizable 
w amounts to a condition on 4>i and wq, and one winds 
up with a one-parameter family of allowed solutions to 
(|B4p . Again I restrict attention to solutions with w > 
everywhere on the expectation that solutions where w 
has nodes are never thermodynamically preferred. 



Vp V 2pi 



q=l 



$1 

Pi 



J 
P 



Wi 

Pi 



(B7) 



I found that pi is never greater than about 3.65, and 
as this value is approached, Wq 0. Using (|B7p . this 
translates into a critical temperature 



Tr « 1.97a 



(B8) 



(See figure [2l) I constructed solutions to (|B4p corre- 
sponding to a range of temperatures ■^Tc < T < Tc- 
A summary of their thermodynamic behaviors is shown 
in figure [3l 

At lower temperatures, the numerical problem be- 
comes more difficult, partly because a larger range of 
r must be integrated over. Also, for smaller values of T, 
extra solutions exist where w has nodes. This is problem- 
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atic because shooting algorithms generally rely on numer- are multiple zeroes, one cannot easily control which of 
ically finding zeroes of a merit function, and when there them one will find. 
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